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Abstract. - A method is proposed for obtaining a systematic expansion of thermodynamic 
functions of spin systems with large spin S in powers of 1/S. It uses the cumulant technique 
and a coherent-state representation of the partition function Z. The expansion of Z in terms 
of cumulants yields an effective classical Hamiltonian with temperature-dependent quantum 
corrections. For the Heisenberg quantum Hamiltonian, they have a non-Heisenberg form. The 
effective Hamiltonian can be solved by methods familiar for classical systems. 

In recent years considerable interest in condensed matter physics has been focusing on the 
properties of new magnetic materials, such as magnetic molecules with a large effective spin 
5, or intcrmctallic compounds containing magnetic layers or chains. For the latter, most of 
the theoretical studies deal with 5=1/2 systems which, in the one-dimensional case, allow 
for analytical solutions using special methods such as the Bethe Ansatz or Jordan- Wigner 
fermion representation of spins. On the other hand, constructing expansions in powers of 1/5 
for systems with large values of S, which are not constrained to special lattice geometries or 
to low temperatures, has not received sufficient attention. 

The aim of this paper is to show that quantum corrections to the classical thermodynamics 
of spin systems with large or moderate 5 can be calculated to an arbitrary order in the 
quasiclassical parameter 1/5. This is achieved by using a cumulant expansion technique 
and working in the basis of spin coherent states. The results which can be obtained by our 
method are valid for temperatures down to the quantum energy scale J/5, where J = J5 2 
is the classical energy scale and J is the exchange interaction. They are complementary to 
those following from the spin-wave theory (SWT) based either on the Holstein-Primakoff or 
path- integral formalisms, which works in the range T < J (see, for example, ref. [El). 

Averages or matrix elements (...) can be expressed through cumulants (. . .) c as follows 

(A) = {A) c , (A l A 2 ) = (A^Atf + (A 1 )(A 2 ), (1) 
(AxA 2 A 3 ) = {A X A 2 A 3 ) C + (A 1 )(A 2 A 3 ) C + (A 2 )(A 1 A 3 ) C + (A 3 )(A 1 A 2 ) C + (A 1 )(A 2 )(A 3 ), 

(*) cnls.lanl.gov/~kladko/; kladko@lanl.gov 

(**) www.mpipks-dresden.mpg.de/~garanin/; garanin@mpipks-dresden.mpg.de 
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etc., where Ai are classical stochastic variables or quantum-mechanical operators. The aver- 
aging above is performed over a classical distribition function or weighed over quantum states. 
A detailed discussion of cumulants, with an emphasis on quantum systems, can be found in 
refs. |3j. Cumulants can be obtained by differentiation of a generating function, i.e., from 



(A,... A N ) C = A • • • -f- ln(e AlAl . . . e x » A ») 
aX% o\n 



(2) 

A 1 = ...=A JV =0 



in contrast to averages (...), which are given by a similar expression without logarithm. By 
setting Ai = A, multiplying by A" /n\, and summing over n — 1, . . . , oo one obtains the identity 

\n(e XA ) = (e XA - 1} C . (3) 

Let us now point out a fundamental relation between cumulants and a quasiclassical theory. 
It is a characteristic property of quantum mechanics, that the expectation value of a product 
of two observables Ai and A2, evaluated with respect to some quantum state is generally 
distinct from a product of the expectation values of A\ and A2, i.e., 

(A 1 A 2 ) c = (A 1 A 2 )-(A 1 )(A 2 )^0, (4) 

whereas equality is achieved in the classical limit. Thus one expects that the theory formulated 
in terms of cumulants will yield a quasiclassical expansion. 

To this end, we rewrite the partition function Z of the system in the basis of spin coherent 
states |n), i.e., maximum-spin states pointing in the direction of the unit vector n. In this 
basis, the cumulants have a simple form, especially if the spin-operator components S z (axis z 
along n) and S± = S x ± iS y are used. First, a cumulant vanishes if it is of the form (. . . S+) e , 
(. . . S+S-) c , (. . . S z ) c , etc. Non-vanishing cumulants are 

(S + Sl l S-) c = 2S(-l) n , 

(S+Sl l S + STS-S r z l 'S-) c = -AS{-l) n+n \-2) m , 

(S + S+S+SS-S-) C = 3(S+S + S-S+S-S-} C = 245, (5) 

etc. The above results can be obtained recurrently, using eq. (|l|) and commutation relations 
[S" z ,5-] = — 5- and =25*2, which remain valid inside cumulants. For scaled spins 

S/S, each nonvanishing cumulant containing n spin operators scales as 1/S n ~ 1 . 
In the overcomplete basis of spin coherent states, the unit operator is resolved as 



25+1 

47T 



dn|n)(n|. (6) 



Thus, for a single-spin system, the trace of an operator A over any complete orthonormal basis 
I to) may be written in terms of spin coherent states as follows |4| 

* 2S +1 f , 2S + 1 f 
tiA = y](m\A\m) = — / dn^(m|A|n)(n|m) = — / dn(n| J 4|n). (7) 

rn m ** 

The partition function for a Hamiltonian H becomes 

Z = trexp(-l3H) = —£— / dn<n| exp(-/3tf)|n>, 0=^ (8) 

The matrix element inside the integral can be written with the help of eq. (^) as 

(n|exp(-/?£0|n) = exp[-/?W(n, /?)], /3H{n, 0) = (n|l - exp(-/3 J ff)|n) c , (9) 
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where W(n, f3) can be considered as an effective quasiclassical Hamiltonian. In order to find 
an expansion of Z in powers of 1/5, we expand the expression for Ti. in a Taylor series: 

H(n,(3) = (n\H\n) c - ^(n\HH\n) c + ^(n\HHH\n) c + ... (10) 

Here the first term on the right-hand side is the classical energy of the spin. As seen from 
eq. (^|), increasing powers of 1/S appear in each order of this expansion. If the quantum 
corrections in eq. ( |10|) are small, then in most cases one can further expand exp(— f37i) in 
powers of 1/S. In this case the statistical weights of different spin orientations are described 
by a purely classical Hamiltonian, whereas quantum effects manifest themselves in corrections 
to the density of states. 

As an example, let us consider the problem of a spin in a magnetic field. 

H = -s-n. (ii) 

For this model, the effective Hamiltonian H(n,/3) can be computed analytically from eq. (||), 
if one expresses the coherent states in the basis of states with definite projections of S on H. 
This yields 



V- -Hn V r^ H (25)! A + COsfl \ S+m / l-CQ S t 

H--Tln^e {s + m)KS _ m)[ [—^—) [—^—) (12) 



m=-S 

where 6 is the polar angle. To obtain the cumulant expansion of Ti, it is convenient to represent 
the spin operator S in cq. ( JTo| ) in the coordinate system with the z axis along the vector n 

S = nS 2 + n+S+ +n_S'_, n± = (n x =p m. y )/2, (13) 

where n x and n y are appropriate transverse basis vectors. Then cumulant averages of different 
spin components are done with the use of eq. . To order 1/S 2 one obtains 

W = -n-h- A[nxh] 2 +I ^(n.h)[nxh] 3 + o(|^, (14) 

where h = 5*11 is the scaled magnetic field. A similar quasiclassical Hamiltonian for a particle 
in a potential V(x) has been recently obtained with the help of cumulants in ref. ||. 

The function Z can be evaluated from eq. @ as for a classical system by doing the integral 
over the directions of the vector n. Physical quantities such as the energy U and magnetization 
m can be obtained by differentiating Z with respect to the appropriate parameter. Direct 
evaluation of observables as classical averages (. . .} must be done with care. For the energy 
one finds 

whereas the scaled magnetization m = (S)/S is given by 

m = = (n >, »=-^="- S [»x[nxh]]+O^J. (16) 

One notices that not only the effective quasiclassical Hamiltonian TL but also the averaged 
quantities contain quantum corrections. 

Next we have to find the limits of applicability of the quasiclassical expansion. It is clear that 
the method will break down at low temperatures. At such temperatures the spin fluctuates by 
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Fig. 1. - Quasiclassical expansions of energy and heat capacity for the spin-in- a- field model, 
eq. (0). 



a small angle from the direction of the field, such that the deviation of the zeroth-order term 
of f¥H [see eq. (|l4])] from the ground state is f3hQ 2 ~ 1. The second- and third-order terms 
of 0TL become (3h6 2 (f3h / S) and f3h9 2 ((3h/S) 2 , respectively. Thus the quasiclassical expansion 
converges for /?/i/5< 1, i.e., for temperatures T>h/S which may be well below the classical 
energy scale h for large spins. For comparison, a usual high-temperature series expansion 
converges for T>h. The quantum energy scale h/S = H is defined by the separation of 
energy levels of the spin in the field. Replacement of discrete energy levels by a continuum of 
classical spin orientations leads to a temperature dependence of the effective Hamiltonian TL or 
of the density of states, if exp(— f3TL) in eq. (||) is expanded in powers of 1/5. This expansion is 
justified under the same condition [3h/S< 1. For (3h/S 3> 1, i.e., T <C h/S, the spin remains 
in its ground state, and a continuous quasiclassical expansion becomes invalid. The form of TL 
in this region is obtained by setting m = 5 in eq. @: TL = -SH - 25Tln[(l + cos0)/2]. 

If exp(— (3TL) is expanded in powers of 1/5, the breakdown of the quasiclassical expansion 
occurs at lower temperatures and the overall results are better. In this case an analytical 
calculation is possible. Using eq. (|l4|), one obtains to order 1/S 2 

bZ a M2S) + b + mi + +0 (Js)- (-) 

where ^ = j3h and B = B(£) = coth £ — 1/£ is the Langevin function. Quasiclassical expansions 
of different order in 1 / S for the energy and heat capacity of the spin-in-a-field model are shown 
in fig. 1. Expansions of orders higher than two have been obtained directly by expanding the 
exact partition function Z = sinh{£[l + 1/(25)]}/ sinh[£/(25)] in powers of 1/5. It is seen that 
already the second quantum approximation nearly quantitatively reproduces the maximum of 
the heat capacity, although the spin value 5 = 5/2 is not very high. 

As a second example, let us consider the model of two ferro- or antiferromagnetically coupled 
spins 

ff = +JSi-S 2 . (18) 

This model is exactly solvable. The eigenstates labelled by the total spin n — 0, 1, . . . , 25 
have the energy E n — +J[n(n + l)/2 — 5(5 + 1)] and degeneracy 2n + 1. The ferro- and 
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Fig. 2. - Quasiclassical expansion for the energy and heat capacity of the ferro- and antiferro- 
magnetic two-spin models, eq. ( |l8| ) . For the antiferromagnetic coupling, the 1/S approximation 
reproduces the exact ground-state energy, Uq = —1.4 J, whereas the 1/S 2 approximation gives 
U = -1.442666 J. 



antiferromagnetic ground states are E2S = —J and Eq = —J(l + 1/S), respectively, where 
J = JS 2 is the classical energy scale. The separation of levels near the ground state is 
AE ~ JS — J/S for the ferro- and AE ~ J = J / S 2 for the antiferromagnet. 

Derivation of the effective Hamiltonian 7i(ni,ri2) for this model, where nx,2 are the 
coherent-state vectors at both sites, proceeds along the same lines. Eqs. (|jH([lO|) retain their 
form, only the matrix elements are defined with respect to both coherent states ni and ri2. The 
latter slightly modifies the routine of calculating the cumulants of the Hamiltonian H : In eq. 
( [l0| ) cumulants of H are rewritten explicitly according to eq. ([!]), then independent averages 
on sites 1 and 2 are taken. After that spins Si and S2 are expressed in their own coherent-state 
coordinate systems and the matrix elements of their components <Si Ql and S2a 2 are computed 
with the help of cumulants. In a cumulant expansion of eq. (]l0|), this modification results in 
the appearance of terms carrying higher powers of 1/S, in addition to the leading terms of 
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order (f3/S) n . To order l/S 2 one obtains 

8 f 2 Rl 2 
H = TJ ni • n 2 - ^[1 - (m ■ n 2 ) 2 ] - - ni . n 2 ] 2 

t S [i - (ni ■ n2)2][i ~ 5ni ■ ns] + (Jp) • (i9) 

Note that quantum corrections have a non-Heisenberg form. As follows from the cumulant 
expansion of cq. (|l0|), an additional pair ni 112 is associated with each further power of /3. Eq. 
( fl9| ) results in the expansion for In Z of the form 

In Z = 2 In(2S) + In ^_ j + + ^ + O j , (20) 

where £ = /3J. 

Quasiclassical expansions for the energy and heat capacity of the two-spin model are shown 
in fig. 2. They are obtained by differentiating eq. (po|). For ferromagnetic coupling the results 
resemble those for the spin-in-a-field model, although here quantum corrections are about two 
times larger. For that model, the leading asymptote U = — /i£(l + l/S)/3 at high temperatures 
is recovered already at order 1 /S, whereas the proper leading asymptote U = — J£(l + 1/ S) 2 /3 
for the two-spin model is obtained only at order 1/5 2 . 

For antiferromagnetic coupling and to order l/S 2 , U and C — dU/dT do not diverge 
as T — > 0, because the leading £ 2 terms in eq. ( |20| ) cancel. The accuracy of the l/S 2 
approximation over a wide range of temperatures is quite impressive. It is tempting to explain 
this by the fact that for antiferromagnetic coupling the separation of the low-lying energy 
levels is by a factor of l/S <C 1 smaller than for ferromagnetic coupling, as was mentioned 
above. Therefore the properties of the system are expected to be "more classical" . 

Furthermore, one could speculate that for many-spin systems without an energy gap the 
quasiclassical expansion works even better. However, it is not true. The quasiclassical 
expansion without a resummation does not reproduce rational powers of T, which appear, 
e.g., in linear spin- wave theory. Therefore, it should break down below the quantum energy 
scale JS = J/S. Note that the SWT based either on the Holstein-Primakoff or path-integral 
formalisms (see, e.g., ref. jj]) works in the range T<^J. Its results should overlap with 
those of the cumulant expansion in the range J/S <T < J. It is interesting to speculate 
how j3/ S corrections of the cumulant expansion go over to 1/5 corrections of the SWT (e.g., 
for antiferromagnets) when T < J/S. One would expect that the leading terms (fl/S) n of eq. 
( fL9| ) convert to the leading order of the SWT, whereas (3/S 2 and other subdominant terms 
result in l/S spin-wave corrections. 

We have shown that for spin systems with large S an application of the cumulant technique 
together with spin coherent states leads to effective, classical-like Hamiltonians with quantum 
corrections taken into account. These Hamiltonians can be further treated by a number of 
methods for classical spin systems, such as exact solutions, Monte-Carlo simulations, or the 
classical-spin diagram technique [pj, which is again based on cumulants. 



REFERENCES 



[1] Auerbach A. , Interacting Electrons and Quantum Magnetism, (Springer, Berlin) 1994. 
[2] Fulde P. , Electron Correlations in Molecules and Solids, (Springer, Berlin) 1995. 
[3] Kladko K. and Fulde P. , Int. Journ. Quant. Chem., 66 (1998) 377; |cond-mat/9709044| Eprint, 
1997. 



K. kladko et al: CUMULANT expansion for systems with large spins 



7 



[4] Lieb E. H., Comm. Math. Phys., 31 (1973) 559. 
[5] GARANIN D. A., Phys. Rev. B, 53 (1996) 11593. 



